Dynamics of spectrally truncated inviscid turbulence 
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The evolution of the turbulent energy spectrum for the inviscid spectrally truncated Euler equa- 
tions is studied by closure calculations. The observed behavior is similar to the one found in Direct 
Numerical Simulations (Cichowlas, Bona'ititi, Debbasch and Brachet [Phys. Rev. Lett. 95: 264502, 
2005]). A Kolmogorov spectral range and an equipartition range are observed simultaneously. Be- 
tween these two ranges a "quasi-dissipative" zone is present in the kinetic energy spectrum. The 
time evolution of the wavenumber that marks the beginning of the equipartition range is analyzed 
and it is shown that spectral nonlocal interactions are governing this evolution. 



Recently renewed interest was shown in the evolu- 
tion of inviscid turbulence. Cichowlas, Bona'ititi, Deb- 
basch and Brachet [l[ (in the following abbreviated as 
CBDB), reported results of Direct Numerical Simulations 
(DNS) of the spectrally truncated 3-D incompressible Eu- 
ler equations, with resolutions of 256 3 , 512 3 , 1024 3 and 
1600 3 wave-modes. In principle, the presence of a high- 
frequency spectral truncation allows the flow to reach a 
thermal equilibrium, i.e. an equipartition of energy over 
wave-vectors. The spherically averaged energy spectrum 
E{K) associated to this equilibrium has a K 2 wavenum- 
ber dependence. This is effectively observed in the CBDB 
DNS, where at high wavenumbers the energy spectrum is 
found to increase following a K 2 dependence. It was also 
observed in CBDB that during the transient towards the 
equilibrium, E{K) followed a K~~ 5 / 3 Kolmogorov scaling 
in the intermediate range. In between these two power 
law regions, a "quasi-dissipative" zone is found in which 
the energy spectrum falls off faster than K~ 5 ^ 3 until a 
local minimum of the spectrum is reached. 

The existence of a K 2 range in the spectrum is not 
new and equipartition was discussed long ago in the lit- 
erature. 0, [I] The observation of a Kolmogorov scaling 
in the energy spectrum of the truncated Euler equations 
is more interesting as it shows that viscous dissipation is 
not mandatory for the build up of an inertial zone but 
that a flux of energy to the small scales, where it accu- 
mulates, is sufficient. It is also an interesting observation 
that between the two coexisting power laws, a dip or 
quasi-dissipative zone, is detected. 

The Eddy-Damped Quasi-Normal Markovian theory 
(EDQNM) H H is known to be compatible with both 
the equipartition of kinetic energy 0, Q and the exis- 
tence of a K~ 5 ^ 3 inertial range. Calculations at higher 
resolution than DNS can be performed at a much lower 
computational cost. EDQNM therefore appears as an 
adequate tool to investigate the high-resolution spectral 
dynamics of the incompressible Euler equations. 

In the present study it is first shown that the spectral 
behavior observed in DNS can be reproduced satisfacto- 
rily by this relatively simple statistical closure. Results 
for resolutions currently unattainable by DNS are ana- 
lyzed. Eventually, it is shown that the dynamics of the 



"quasi-dissipative" zone are governed by highly nonlocal 
triad interactions. 

In the case of the EDQNM closure, the evolution equa- 
tion for the energy spectrum corresponding to the Euler 
equations is the Lin equation without viscosity: 



dE(K,t) 
dt 



= T NL (K,t) 



(1) 



in which the nonlinear transfer TVl is expressed as: 

T NL {K,t) = J J Q KPQ (xy + z 3 ) [K 2 PE(P,t)E(Q,t) 

dPdQ 



-P 3 E(Q,t)E(K,t)] 



PQ 



(2) 



In equation @ , A is a band in P, Q-space so that the 
three wave- vectors K,P,Q form a triangle. x,y, z are 
the cosines of the angles opposite to K, P, Q in this tri- 
angle. The characteristic time Qkpq is defined as: 



e 



KPQ 



1 - exp(-(r) K + VP + Vq) x *) 



T)K + VP + VQ 

in which 77 is the eddy damping, expressed as 



Vk = A 




S 2 E(S, t)dS. 



(3) 



(4) 



For A the classical value 0.36 is retained corresponding 
to a value for the Kolmogorov constant of 1.4. The 
spatial resolution for the computations reported here is 
14.2 wavenumbers per octave. This significantly high 
wavenumber density was selected in order to satisfacto- 
rily capture the nonlocal interactions that will be shown 
to play an important role in the following. The effect 
of the truncation of the domain at a cut-off wavenumber 
Kmax is introduced by omitting in (|2|) all interactions 
involving wavenumbers larger than K max . 

In the present paper we analyze a freely evolving invis- 
cid velocity field (the term decaying would be misleading 
since no energy is dissipated). The initial spectral en- 
ergy distribution is localized at small wavenumbers. The 
initial energy spectrum is: 



E(K,0) = BK s e- 2K2 / K * 



(5) 
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with s — 4 and B determined so that the total kinetic en- 
ergy is equal to unity. In this spectrum Kl is a wavenum- 
ber characterizing the initial integral length scale, and 
the computations are run on a spectral domain extending 
from a minimum wavenumber Kq to the high frequency 
cut-off K max . These initial conditions arc different from 
the ones used by CDBD, who started from a Taylor Green 
vortex, so that the short time evolutions may be differ- 
ent in the two studies. A first set of computations was 
performed with K = K^. In this case the maximum of 
the energy spectrum is situated at the first wavenumber 
and the K s part of the spectrum is absent. This set has 
the advantage of permitting comparisons with the DNS 
of CBDB in which the energy maximum is also located at 
the low wavenumbers at the beginning of the computa- 
tions. Calculations are performed with ratios of smallest 
to largest wavenumbers K max /Ko = 85, 171, 341, 533. 
These resolutions are equivalent to the 256 3 , 512 3 , 1024 3 
and 1600 3 DNS calculations of CBDB respectively. Cal- 
culations were also performed at higher resolutions, that 
would correspond to DNS on 4096 3 , 8192 3 , 16384 3 and 
32768 3 grids. In the following the different calculations 
will be denoted by 256, 512, etc. Finally we report the 
results of a calculation starting with a spectrum having a 
K A low wavenumber zone. In this case Kl/Kq = 10 and 
K max /Ko — 10 4 . In the following the time is normalized 
by an initial eddy turnover time defined by (MKl) , in 

which U = (§/ E{K)dK) 1/2 . 

The results for the time evolution of the 1600 calcula- 
tion are shown in figure 1. Also shown are the spectra at 
t = 2.4 for different resolutions. Comparison with Figure 
1 of reference [l[ shows that the present closure results 
are very similar to the DNS results. The results of the 
calculation starting from an initial spectrum with a K A 
low wavenumber range are shown in figure 2. 

In all calculations the nonlinear interactions create a 
^--5/3 mer tj a i range. The existence of the inertial range 
is particularly clear in the results at high resolutions. In 
this range the energy spectrum scales as: 



E(K,t) = C^e(t) 2/3 ^" 5/3 , 



(6) 



in which e(t) should not be interpreted as a molecular dis- 
sipation which would be zero in the present case, but as 
an energy flux in the spectral cascade. This range elon- 
gates with time towards higher wavenumbers until the 
end of the inertial range meets the maximum wavenum- 
ber of the truncated domain K max . Then a thermalized 
range, or equipartition spectrum, 

E(K,t) = A(t)K 2 , 



(J) 



starts to build up at large wavenumbers. This thermal- 
ized range later extends from large to smaller wavenum- 
bers, and A(t) is an increasing function of time. The 
thermalized energy eth associated with this range is also 
increasing with time. By integrating the equipartition 
range it is found that A(t) is related to eth by the expres- 
sion A(t) w 2>ethi 'Kmax- For long times the thermalized 
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FIG. 1: Top: time evolution of the energy spectrum E(K) for 
calculation 1600. The solid lines indicate the K~ 5 ' 3 and K 2 
spectral slopes. Middle: spectra for the calculations 256, 512, 
1024 and 1600 at t = 2.4. Bottom: spectra for the calculations 
4096, 8192, 16384 and 32768 at t = 2.4. 
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FIG. 2: Time evolution of an energy spectrum with a K 4 low 
wavenumber range. 
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energy approaches the (invariant) total energy and A(t) 
tends to a constant value. 

The coexistence of a Kolmogorov range and of an 
equipartition zone was already observed by Connaughton 
and Nazarenko p] using the diffusion approximation 
of turbulence proposed by Leith [8j. The "quasi- 
dissipative" zone is however absent in their results. In 
agreement with the DNS results of CBDB, the "dissipa- 
tive" range is present in the EDQNM results. Leith's 
model was shown to be a useful and easily maneagable 
model to obtain understanding of various problems re- 
lated to turbulence (e.g. scalar mixing, two-dimensional 
turbulence Its character is however spectrally 

local. EDQNM, on the contrary, basically relies on the 
triadic interaction between wave-modes and is therefore 
by essence nonlocal. The fact that the " dissipative" range 
is reproduced by the EDQNM closure, while it is absent 
in the results obtained with Leith's diffusion model 0, 
indicates that the nonlocal interactions between the ther- 
malized modes and the modes in desequilibrium might be 
responsible for the existence of this zone. This issue will 
now be addressed. 

As in CBDB we introduce K t h{t) as the wavenum- 
ber corresponding to the minimum of the energy spec- 
trum. Kth(t) characterizes the beginning of the absolute 
equilibrium zone. From the equations © and |7|). it is 
straightforward to calculate the intersection of the two 
zones, yielding a first estimate for K t h(t): 



K p i{t) 



e(t) 



2/11 



#9/11 

max i 



(8) 



The subscript pi indicates that a K~ 5 / 3 power law iner- 
tial range is presumed reaching the equipartition range 
at K — K t h(t). This estimate was first proposed in Ci- 
chowlas el ai.pj]]. Another possible estimate is: 



K d (t) 



e(t) 



1/4 



#3/4 



(9) 



It was also proposed in reference [Tl| (see also CBDB). 
The physical argument behind this estimation is the ex- 
istence of a "quasi-dissipative" region in the spectrum, 
occurring between the Kolmogorov and equipartition 
ranges. A simple way to find Kd by dimensional argu- 
ments is to postulate that Kd is the inverse of a Kol- 
mogorov scale 



1 



i/4 



(10) 



built on the energy flux e and an eddy viscosity v t . As- 
suming v t determined by the thermalized energy and the 
wavenumber bound of the spectral domain K max , 



Vl 



K„ 



(11) 



§§§ is immediately found. Introducing K~^ ax as the 
lengthscale on which to build the effective viscosity v t 
is an assumption typically nonlocal in wave-space. 
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FIG. 3: Evolution of the parameter K t h(t)- Top: K t h com- 
pared to K p i (equation (JSj ) . Bottom: K t h compared to Kd 
(equation ([9])). 



The nonlocal nature of estimate ((9|) can be enlightened 
using the EDQNM equations and performing a nonlocal 
expansion a la Lesieur and Schertzer [l2| . Developing the 
nonlinear transfer ([2]) with respect to Kj Ki and assuming 
K « Ki and E{K) < E(Ki) yields: 
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:K 4 







KPP- 



T NL {K) = -2u t K 2 E(K) 
E(Pf 



K, 



P 2 



-dP + <D[KE(K)} 3/2 (12) 



Focusing on the first term that has the form of an eddy 
viscous term with: 



v t = T7 



K„ 



Ki 







KPP 



5E(P) + P 



dE{P) 
dP 



dP 



(13) 



and assuming E(Ki) in the equipartition range, it is read- 
ily found that v t is given by expression (|11[) . The fact 
that expression (flTj) . and therefore estimate (9), can be 
found by a non local expansion of the tranfer term, un- 
derlines the non local nature of the physical mechanisms 
associated with this estimate. It has to be noted that ex- 
pression ^ was proposed by CBDB following a different 
argument. 

While in DNS it is not easy to decide which estimator 
for K t h is the more appropriate, the relative cheapness 
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FIG. 4: Evolution of K t h{i), compared to the two estimates 
K p i(t) and Kd(t) (equations ([8} and J9} respectively) for the 
K 4 case. 

of EDQNM calculations allows for high resolution calcu- 
lations that help to choose between the two estimates. 
In figure 3 the values of K t h found with the closure are 
compared to the estimates © and (0). The dissipative 
estimate is clearly found to be the more appropiate since 
it leads to a better collapse of the results on the different 
runs as well as to a more pronounced horizontal plateau 
of the plotted ratios. In figure 4 it is confirmed that also 
in the K A case the dissipative estimate of K t h{t) is supe- 
rior to the power-law estimate. It is interesting to point 
out here that the opposite conclusion would be obtained 



with a model that excludes the nonlocal interactions as 
in Connaughton and Nazarenko pj. 

The main results of the study can be summarized as 
follows. The EDQNM closure reproduces the behavior 
observed in DNS, i.e. an energy spectrum containing a 
Kolmogorov inertial range, a dissipation range and an 
equipartition range. The dissipation range was shown to 
be created by nonlinear interactions with the modes in 
equipartition. An effective eddy viscosity can be defined, 
based on the most energetic modes in the equipartition 
zone and the cut-off wavenumber. The non local charac- 
ter of this eddy viscosity was verified by expanding the 
non-local contributions to the nonlinear transfer in the 
EDQNM formulation as a function of the wavenumber 
ratio and retaining the leading order term. 

Although the problem of truncated inviscid turbulence 
remains somewhat artificial and far from real world ap- 
plications, some of the spectral behaviors studied here 
may help understanding mechanisms that are connected 
to Navier-Stokes turbulence. It is for example interest- 
ing to note, in the framework of Large- Eddy Simulation, 
that the turbulent cascade of energy at large scales is not 
significantly altered by the fact that energy is not, or not 
properly, dissipated. 

We acknowledge stimulating interaction with Marc- 
Etienne Brachet in the framework of the GDR "Structure 
de la Turbulence et Melange" of the Centre National de 
la Recherche Scientifique (CNRS), which provided the 
starting point for the present work. 
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